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We offer theoretical explanations for some recent observations in numerical simulations of quantum 
random walks (QRW). Specifically, in the case of a QRW on the line with one particle (walker) and 
two entangled coins, we explain the phenomenon, called "localization", whereby the probability 
distribution of the walker's position is seen to exhibit a persistent major "spike" (or "peak") at 
the initial position and two other minor spikes which drift to infinity in either direction. Another 
interesting finding in connection with QRW's of this sort pertains to the limiting behavior of the 
position probability distribution. It is seen that the probability of finding the walker at any given 
location becomes eventually stationary and non- vanishing. We explain these observations in terms 
of the degeneration of some eigenvalue of the time evolution operator U{k). An explicit general 
formula is derived for the limiting probability, from which we deduce the limiting value of the height 
of the observed spike at the origin. We show that the limiting probability decreases quadratically for 
large values of the position x. We locate the two minor spikes and demonstrate that their positions 
are determined by the phases of non-degenerated eigenvalues of U(k). Finally, for fixed time t 
sufficiently large, we examine the dependence on t of the probability of finding a particle at a given 
location x. 
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PACS numbers: 03.67.Lx, 05.30.-d, 05.40.-a, 89.70.+C 



I. INTRODUCTION 

Discrete-time random walks employing a quantum coin 
originated in 1993 with the publication of [jy]. Since then, 
inspired by the promise of applications to the develop- 
ment of super-fast algorithms, especially after the publi- 
cation of J:, 4, i5j] in the years 2000 and 2001, this field of 
research has flourished immensely. 

A natural spatial setting for the evolution of a discrete- 
time random walk is an infinite linear lattice. In the sim- 
plest classical example, the walker starts at the origin 
and subsequently, depending on the outcome (heads or 
tails) of the toss of a fair coin, strolls one step (spatial 
unit) either to the right or to the left. In the quantum 
context, the tossing of the coin is modeled by the action 
on the overall state function of a "coin operator" while 
the movement of the walker is modeled by the action of 
a "shift operator". The observable aspects of the sys- 
tem are captured by the eigenvalues of these operators. 
The characteristically quantum-mechanical phenomenon 
of superposition of states gives rise to spatial probability 
distributions quite unlike those observed in classical ran- 
dom walks. Typically, sharp spikes of high probability 
density are observed at specific locations and troughs of 
near zero probability elsewhere. 

A number of cases have been studied in some depth. 
For instance, the case of one particle and one coin is 
treated in p, H, [sl. [ill. [l4| : two entangled particles and one 
coin in [l7[ ; one partilce and many coins in [9|, ll9| , and 



one particle and two entangled coins in [l5| , etc. Very re- 
cently, Venegas-Andraca et al. [l5| showed, by way of nu- 
merical simulations, that the walker tends to persist with 
high probability at the initial position as evidenced by a 
very prominent "peak" or "spike" at the origin. The oc- 
currence of spikes and similar phenomena (called localiza- 
tion) also has been reported in two-dimensional Grover 
walks [7] . llOl I12I |20| , in multi-state quantum walk on a 
circle [l3J, in one-dimensional quantum walks driven by 
many coin s [p and in one-dimensional three-state quan- 
tum walk [lyl • 

In this paper, we concentrate on the case of a discrete- 
time QRW on a linear lattice with two entangled coins, as 
proposed by Venegas-Andraca et al. [15^. In this frame- 
work, abbreviated "2cQRW" , the controlling behavior of 
the two entangled coins is modeled by a certain 4x4 
matrix which is defined as the second-order tensor-power 
of a 2 X 2 unitary matrix. As in [1|| , we show that the 
occurrence of localized spikes at the origin and elsewhere 
refiects the degeneracy of some eigenvalue of the time 
evolution operator U{k). An explicit formula for the lim- 
iting spatial probability density enables us to specify the 
height of the spike at the origin. For sufficiently large val- 
ues of the position x, we find that the limiting probability 
density decreases quadratically. This is in sharp contrast 
to the result obtained by Inui et al for a three-state QRW 
[l6|, where it was found that the limiting spatial proba- 
bility density decreases exponentially for large x. 
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II. FORMULATION OF QUANTUM RANDOM 

WALKS WITH ENTANGLED COINS ON THE 

LINE 

A. Quantum random walks on the line with two 
entangled coins 

We proceed to define the elements, as formulated in 
[15J . of a random walk on a discrete linear lattice gov- 
erned by a pair of entangled qubits. The definitions and 
corresponding notations are analogous to those for the 
single-coin (one-qubit) framework outlined in [2l|. For 
brevity of exposition, as in [l5j , the two-coin, or more pre- 
cisely, two-qubit framework is abbreviated by the name 
"2cQRW". 

Akin to the one-qubit case, the coin space of the 
2cQRW framework is the Hilbert space Hec ( "ec" for en- 
tangled coin) spanned by the orthonormal basis {\j)',j G 
Be} where Be = {00,01,10,11}. The position space is 
the Hilbert space Hp spanned by the orthonormal basis 
{|a;);a:; S Z}. The "overall" state space of the system is 
H = Tip (8) Tiec, in terms of which a general state of the 
system may be expressed by the formula: 



V' = X! X! V'(2^,j)|a;)«' Ij). 



xeZjGBa 

The spatio-temporal progression of the 2cQRW is gov- 
erned by an "evolution operator" U, which is composed 
of a coin operator A^c and a shift operator S. 

In the 2cQRW context, the coin operator is defined as 
the tensor product of two single-qubit operators 

Aec^ A® A, 

where the unitary operator A acts on a single-coin space 
(see P). 
As in [1^1, the shift operator is given by: 



s = mm ® ^ N + i){i\ + ioi)(oi| ® Yl i^)(*i 

i i 

-f iio>(ioi ® Y. N)(*i + iii)(iii ® E I* - i)(*i- w 

i i 

The journey of the particle (a.k.a. walker) along the 
line is driven by a stochastic sequence of iterations of S. 
At every time step of the walk, depending on the state of 
the coin, the walker strolls either one spatial unit to the 
right or to the left or stalls at that step. Note that the 
walker moves only when both coins reside in the |00) or 
1 11) state. Otherwise, the walker stalls at that step. 

Let / denote the identity operator on Tip. Then, in 
terms of A^c and S*, the total evolution operator U is 
given by 



U = S{I®Aec). 



Given tpo G H, where 



1, the expression tpt 



t. The corresponding random walk with initial state f/'o 
is represented by the sequence {'0t}o°- 

Let X denote the position operator on Hp, defined by 
X\x) = x\x) and let Vt = E^GzEjeSc Mx,j)\x) ® \j) 
be the wave function for the particle at time t. Then the 
probability pt (x) of finding the particle at the position x 
at time t is given by the standard formula 



Ptix) 






where | • | indicates the modulus of a complex number. 
At each instant t, the eigenvalues of the operator Xt = 
[/''' XU* equate to the possible values of the particle's 
position with corresponding probability pt{x). 



B. Fourier transform formulation of the w^ave 
function for 2cQRW 

As in the single-coin setting treated in [ll|, Fourier 
transform methods can be applied in the 2cQRW set- 
ting to obtain a useful formulation of the wave func- 
tion. Let ^rW = [Mx,i),Mx,2),Mx,^),MxA)f 

represent the amplitude of the wave function, whose 
four components at position x and time t correspond 
respectively to the coin states 00, 01, 10, and 11. 
As usual, the superscript T denotes the transpose 
operator. Assuming that the 2cQRW is launched 
from the origin, then the initial quantum state of the 
system is reflected by the components of ^§'^(0) = 
[V'o(0,l),^o(0,2),V'o(0,3),V'o(0,4)]^ = [ai,a2,a3,aif, 

where J2j=i l"iP = 1- 

To begin, the spatial Fourier transform of ^('^(x) is 
defined by 



^(fc) = ^*r(2^)e' 



xez 



For instance, under this transformation, the initial am- 
plitude is related to its Fourier dual by the formula: 



%%k) = M/-(0). 



(2) 



In general, the Fourier dual of the overall state space 
of the 2cQRW system is the Hilbert space L^(K) (g) Hec, 
consisting of C*^- valued functions: 



0(fc) 



Mk) 
Mk) 
Mk) 



(3) 



subject to the finiteness condition 



^illi^ 



^211^2 



hWh 



»4||l2 



< oo. 



Thus, given the initial state '^^'^{k), the Fourier dual of 
the wave function of the 2cQRW system is expressed by 



U i/jQ is called the wave function for the particle at time 



^rik) = C/ec(A:)**-(fc), 



(4) 



where the total evolution operator Uec{k) on L^(K)(; 
is given by 

ik 



Uec{k) = 

Note that Uedk) = 
U{k/2) 











-ik 



A, 



C/(fc/2)(8)C/(fc/2), where 







-ik/2 



A. 



me 



(5) 



(6) 



To wrap up this introductory section, we collect some 
basic facts about the eigenvalues and eigenvectors of the 
operators discussed above. 

Since the matrix A is unitary, we may assume, with- 
out loss of generality, that its determinant is \A\ = e*^, 
where is a real constant. Similarly, for the unitary 
matrix [/(fc/2), we may assume that its eigenvalues are 
Ai(fc) = e*''('=) and AaCfc) = e*(^-''('=», where ry is a 
real- valued differentiablc function of k. Let vi{k) and 
V2{k) denote the corresponding unit eigenvectors. Since 
Uec{k) = C/(fc/2)(8)[/(fc/2), it follows that the eigenvalues 
of Uec{k) are: 



' Ai(fc) = [Ai(fc)]2 = e^v(k) ^here Lp{k) 
K2(k)^\dk)\2{k)^\A\=e'' 
K^(k)^\dk)X2{k)^\A\ 

[ A4(fc) = [\2{k)f 



2r]{k) 



AB 



^i(20-^(k)) 



Correspondingly, the unit eigenvectors are: 

f Viik) = vi{k) <g) vi{k) 
V2{k)=vi{k)(S)V2{k) 
Vslk) ^ V2ik) (g) vilk) 
{ Vilk) =V2{k)®V2{k). 
Finally, in terms of eigenvalues and eigenvectors, the 
wave function '^'^'^{k) may be expanded as follows: 

^(fc) = ui{k)n%^) 

+e"'{V2{k),^%k))V2{k) 
+e'''{V^{k),%'^{k))V^{k) 
+e''{'i'>-v''^'>\Vi{k),^dk))Vdk). (7) 



III. 



SPATIAL PROBABILITY DISTRIBUTION 
OF 2CQRW 



For simplicity of notation, occasionally the explicit de 
pendency on the parameter k of the quantities {V^j^^i 

'^f^ and '^f^ will be suppressed. By inverse Fourier trans 
formation, the amplitude of the wave function of the par 
tide at the position x and the time t is given by 



10 

h 



g-«xfc^e 



h+h 



27r 



(8) 



where 



^1 = /o e-"'=e"'^W(V^i,*g=)Fi II 



J r'^''^ „ — ixkAt8 l\r iT/ec\T/ dk 

h-Jo e e (1/2, Wq )\/2 2;; 



I3 
h 



/c 



27r 



-ixk AtO 



dk 
dk 



(^3,*S^)V^3f 



^^^ ^-^xk^^t{2Q~^{k))^y^^^e^c^^y^ 



dk 
2ix- 



In view of [I3], the asymptotic behavior of the four 
Fourier integrals representing the amplitude components 
/i, /2, /a and /4 of ^('^(x) can be described by the follow- 
ing lemma: 

Lemma 1. Let ^(fc) be an A'^-fold continuously dif- 
ferentiablc function in the interval < fc < 27r. Let 

and Bn(x) = J2n=o^"~^9^''Hb){-x)~'^-\ Then, as 
|a;| — > cxD, we have 



2Tr 



e-''"'g{k)dk = Bn{x) - An{x) + o{x-'^). (9) 



The following theorem justifies the prominent spike at 
the origin observed by Venegas- Andraca et al in [1^ , as 
well as other interesting phenomena. 

Theorem 1. l^et pt{x)^\\^r{x)f =Y.^^^^\i,t{x,j)\\ 
where |].|| denotes vector norm. Letp(a;) — Yivat^ryo Pt{x) ■ 
If the 2cQRW is launched from the origin with initial 
state ^Q^ and governed by evolution operator Uedk), 
then the limiting probability of finding the walker at |a;) 
is 



p{x) = \\ I {e 



27r 



— ixk \ "^ 



dk. 



^{V,(k),n%k))V,{k)}^f. 

J=2 



In particular 

(i) when x = 0, we have: 



p(0) = 



277 3 



-— ~ rlh 

{}_J(V,-(fc),*-(fc))F,(fc)}-| 

J=2 



(10) 



(ii) as \x\ -^ 00, we have: 



p(^) - illEK^^(o)'*3'(o))y,(o) 



j=2 



{V,i27T),%d2n))V,i2iT)]r + Oix- 



ai) 



Proof. See Appendix A. 

The height of the spike at the origin is quantified pre- 
cisely by Eq. pi]]) . For instance, if we choose, as the ini- 
tial coin state, the BeU state |$+) = -^ (|00) -h |11)), 
then the probability of finding the particle at the origin 
converges to p(0) = 3 — 2-\/2 w 0.171573 (justification 
below). This agrees quite accurately with the graphical 



representation in FIG. 1 (see below) and the observations 
reported in J^. No such spikes are evident in QRWs 
driven by single coins. 

The formula Eq. (jlip for p{x) implies that the limiting 
probability of finding the particle at any fixed location 
does not vanish. This is in sharp contrast to the behavior 
of single-coin QRWs [3, 0, [^, [lJ| , for which the analogous 
limiting probability converges everywhere to zero. 

In [16|], which treats the case of a one-dimensional 
three-state quantum walk, there too the behavior of p{x) 
is observed to be everywhere non- vanishing. However, 
one important difference distinguishes that case from the 
present case. Whereas in [l^ the rate of decay of p{x) 
as |a;| ^ oo is exponential, in the present case the rate of 
decay oi p{x) as |x| ^ oo is quadratic. 

By the proof of Theorem 1 (see Appendix A) , the ev- 
erywhere non-vanishing behavior of p{x) as well as the 
occurrence of the spike at the origin both are due to the 
independence on k of the degenerate eigenvalues of the 
evolution operator Uec- 

We remark that the sequence {p{x)}'!^^q of hmiting 
probabilities does not itself amount to a probability dis- 
tribution. In fact the sum X]x6zP(^) usually is less than 
unity. The following theorem evaluates this sum. 

Theorem 2. If the 2cQRW is launched from the origin 
with initial state ^g'^ and governed by evolution operator 
Uec{k), then 






271 3 

E 

J=2 



\{V,{k),^l^{k))\ 



, dk 

2^' 



(12) 



Proof. See Appendix B. 

To illustrate the above theorem, suppose the 2cQRW 
is launched from the origin with initial state |<I>+) — 
4j (|00) -I- |11)), governed by the coin operator i/®^, 
where H is the 2x2 Hadamard operator. Then, by P^ . 

X]xGZ-P(^) — n/2— 1- This example is discussed further 
in the sequel. 

The third and final theorem of this article provides 
detailed information on the local behavior oi pt{x), the 
probability of finding the particle at time t and position 
X. In preparation for this theorem, we digress for a mo- 
ment to review some prerequisites. 

As in Eq.®, let f{k) denote the phase of the princi- 
pal eigenvalue Ai = ^^viM) associated with the evolution 
operator Uec{k). To ensure sufficient control over the an- 
alytic properties of ip, our choice of the coin operator A^c 
must be restricted accordingly. Specifically, (f must pos- 
sesses no stationary points (see [3|, |l8|) of order greater 
than 1. Moreover, the maximum and minimum values of 
the function p'{k) over [0, 27r] must coincide respectively 
with |(/3'(fco)| and — [(^'(/cq)!, where /co is a zero of f"{k) 
of multiplicity one. 

Provided that f3 is neither nor ^ , the aforementioned 
assumptions are valid at least for the important case of 
the coin operator A(/3)®2 treated in [2l|. If /3 is or f, 
then the resulting QRW is trivial. 



In the special case when f3 = j, our coin operator 
^(f)®^ coincides with the 4x4 Hadamard-based matrix 
jj02^ This example is fully elaborated below. 

Theorem 3. Suppose the 2cQRW is launched from 
the origin with initial state ^g°(fc) and governed by 
evolution operator Uec{k) based on the coin operator 
A^c = A(/3)®^. Let i^(fc) denote the phase function of 
the principal non-degenerate eigenvalue Ai = e*'^^'^-' of 
Uec{k) and let M = \ip'{ko)\, where fco is a zero of 'p"{k). 



(i) When a; = 0, we have 

pt{0) = p(0)+O(t-^). 



(13) 



(ii) If X is an integer within a small fixed distance b of 
±tM then pt(x) = Oitrl). 

(iii) For a fixed e > 0, if t > a: > t{M + e) or t>-x> 
t{M + e), then ptix) = 0{t-^). 

(iv) For a fixed e > 0, if i^ < x < t(M - e) or t^ < 
-X < t{M - e), then pt{x) = 0{t-^). 

(v) If |a;| < t^, and |a;| « 0, then 

Pt{x) = p{x) + 0{t-^). (14) 

(vi) If \x\ < t^, and |a::| w t^ , then pt{x) = 0{t^^). 

According to this theorem, if t is fixed and sufficiently 
large, the distribution oi pt{x) is characterized by a spike 
at the origin and two minor spikes located at ±ii^'(fco). 
As t increases, the two minor spikes shrink in height and 
drift off to infinity in both directions (see FIG. 1), while 
the spike at the origin persists and settles to a specific 
height given by Ea. (|10p . These inferences conffim the 
observations based on numerical simulations reported in 

m- 

Example. For the duration of this example, let us 
agree to adopt the convection whereby any row vector, 
when enclosed in round parentheses, is identified with the 
corresponding column vector enclosed in square brackets. 
For instance: 



(ai,a2,a3) 



0.2 

as 



Suppose the 2cQRW is launched from the origin with 
initial state 1$+) = -i= (|00) + |11)) and governed by the 
coin operator 



TjOiz _ 



-^(|0)(0| + |0)(1| + |1)(0|-|1)(1|) 



which equals the operator A^c — A{f3)'^'^ in Theorem 3 
when (3 = J. 



The total evolution operator is given by Uec{k) 
t7(fc/2)(g)C/(fc/2), where 



uim = ^_ 



ik/2 



ik/2 



g-ifc/2 _g-ifc/2 



Let(p(/fc)=2sin-i(2M|H) 



(15) 
Then the two eigenval- 



ues of U{k/2) are 

Ai(fc) = e"^W/2 

A2(fc) - -e-*'^('=)/2. 

and the four eigenvalues of Uec{k) = U{k/2) U{k/2) 
' Ai(fc) = Ai(fc)2=e*'^('=) 



are: 



A2(fc) = Ai(fc)A2(fc) = -l 
A3(fc) = Ai(fc)A2(fc) = -l 



[ A4(A;) = A2(fc)2=e 
To apply Theorem 3, observe that 
cos(A:/2) 



V'{k) = 
V"{k) 



2-sin2(fc/2) 
sin(fc/2) 



2[2-sin2(A:/2)]! 



so that the zeros of (,c"(fc), both of multiplicity one, occur 

at and 27r. Thus, the maximum value M — H^ oi (p'{k) 

is attained at fc = while its minimum value — -^ is 
attained at fc = 27r. 

For economy of notation, let 



and 



71 (fc) 


= - 


cos 

COS 


1 + ^/1+ C0S2 


k 

2 


72 (fc) 


l-^l + cos^ 


fe 
2 


Ni 


(fc) 


= 2 


- 27i(fc)cos| 




N2 


(k) 


= 2 


-272(fc)cos|, 





in terms of which, the two unit eigenvectors of U{k/2) 
corresponding to the eigenvalues Ai(fc) and X2{k) are: 

t,i(fc) = ^(e''=/2,7i(fc)) 
«2(fc) = ;7Wj(e^'^''72(fc)), 

and the four unit eigenvectors of Uec{k) are: 

f ^i(fc) = WT (e*^e*'=/27l(fc),e^'/'7l(fc),7l(fc)') 
^2(fc) = ^^ (e*^e*'=/272(fc),e^'=/27l(fc),-l) 

^3(fc) = VWm (e*',e*'/'7i(fc),e^'=/2^2(fc),-l) 
^4(fc) = ^ (e*^e*'=/272(fc),e^'=/272(fc),72(fc)2) . 



By Eq.©, with *3=(fc) 



^,0,0,^ 



the wave 



'0 y^J ~ \y'2 ' "' "' 2 

function ^t'^(fc) = C/|g(fc)^g^(fc) may be expressed as a 
sum of four components 



"iZ'ik) ^ Hi{k,t) + H2ik,t) + H3{k,t) + Hiik,t), 

where: 

f i7i(fc, t) = ^2|^ (e-'^ + 7i(fc)2) Vi(fc) 

-^^■-l)T/2(fc) 

i^4(fc, t) = ^^i;;^ (e-^ + 72(fc)2) V^4(fc). 

Applying inverse Fourier transformation, as in Eq.([i 
the four summands of 

*r(a;)=/i+/2+/3 + /4 






(16) 



( T - f^'' 
^1 -Jo 



xk v2' 



|^(e-^'= + 7i(fc)^)Vi(fc)S 



-^2 = Jo e 



c/c v2e 



2vWr?V; 



_ .2, 

-'a -Jo 



^ — ixk v2e 






-ifc 



— ik 



l)V2{k) 



dk 
2-K 



/c 



2ir 



2\/]Vl 

cfc V2e'*'^<''' 
2Af2 



l)V^3(A:)f 



(e-'^-+72(fc)^)X4(fc)f 



To determine the height of the spike at the origin, as 
given by Eq. (fTO]) , we must evaluate the sum of the second 
and third terms of P^ at a; = 0. After some algebraic 
manipulations, using the identity N1N2 — 4 + 4 cos 
and the formulas 



2 fc 
2 



2ir 



1 dk V2 



1 + cos2 I 27r 2 ' 

cosfc dfc 3a/2 



1 + cos2 I 27r 



27r 



sin fc dfc 



0, 



/o 1 + cos2 I 27r 
we arrive at the formula: 
I2 +h = l ((-1)*(2 - 72), 0, 0, (-1)*(2 - V2)) . (17) 

Therefore, the probability of finding the particle at the 
origin converges to 

p(0) = II/2 + /3II 
= 3-2V2 
« 0.171573. 

Our computer simulation (see FIG. 1) accords well 
with this theoretical value as does the simulation result 
0.171242 reported in [15;]. 
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FIG. 1: The position probability distribution pt{x) for a 
2cQRW with initial coin state |$+) = -^ (|00> + |11>) after 
t = 400 steps. 



where 

1 2\/2-2 _ 
Oi ~ = — e 

27r(4-2V2)2 



/ 96 \^r(i) 

\tV2j 3 



1 -^V^ ^ „-ic2TT „-itV2TT^ 

^2 = ::: — : ;=— re e 



96 YT{^) 

tV2y 



27r(4-2N/2)2" " Vi%/2y 3 

Finally, after some further manipulations, we obtain: 

Pt{x) ~ II/1+/4IP 

(6V2)lr(i^2 



e '6 



-ic27r -it\/27r+if || 2 



187r2 



i"3 



(6V2)ir(i)2_ 



67r2 



i"3. 



We proceed to elucidate the destinies of the two minor 
spikes, which plainly are visible in FIG. 1. By Theo- 
rem 3, with M — ^, each of the spikes, respectively, 
is predicted to lurk within a small neighborhood of the 

positions x = ±-^^- It suffices to consider only the spike 

at X = -^t, since the other one may be treated similarly. 
Also, to serve as our small neighborhood, it suffices to 
adopt the interval Jt = [^t - 1, ^t + 1]. 

Let X = ^t + c e Jt- Substituting in each of the 
summands of Eq. p6p this value of x and performing some 
simple algebraic manipulations, we obtain: 



^1 = /o 



2^^,tMk)-^k]^-^k5^ (^-.k ^ ^^(^)2^ y^(^) II 



27Vi 



2^ ^iirt p — ixk 



^2 = /o e 



V2 



2VN1N2 
V2 



(e"''^--l)y2(fc)f 



~ik 



JO 



l)^3(fc) 



dk 
2-K 



e 

4 = /-e^*[-^(fe)-^fe]e-'=^^(e-''=+72(fc)2)l^4(fc) §. 



By Lemma 1, both I2 and I3 decay as 0{t~^). Hence, 
as i ^ 00, their contributions to the value of pt{x) be- 
come eventually negligible relative to /i and /4, whose 
values, as shown below, decay no faster than 0{t^^). 

Without much effort, one verifies that the modified 
phase function ip{k) — ^k possesses a stationary point of 

order 2 at fc = and, by similar reasoning, —ip{k) — ^k 
possesses a stationary point of order 2 at fc = 27r. Thus, 
by the method of "stationary phase" [3, |T^ , the following 
values are obtained for the leading terms of h and I4: 



h 

h 



Ci • (l,V2-l,%/2- l,3-2%/2) 
C2-(l,\/2-l,V2-l,3-2V2) 



IV. RELATED IDEAS FOR FURTHER 
RESEARCH 

As above, let X denote the position operator on the 
position space Tip and let Xt = V^ XU*, where U de- 
notes the evolution operator. A natural and potentially 
more useful statistical description of the evolution of the 
2cQRW as time t ^ 00 can be developed in terms of the 
"normalized" operator \Xi. 

Suppose the 2cQRW, controlled by the coin oper- 
ator iJ®2^ is launched from the origin in the initial 
state *§'=(0) =ai|00) -Ha2|01) -HaallO) +ai\ll), where 
X]i=i I'^jP = 1- -'^°'' y ^ ["Ijl]: Ist 5^{y) denote the 
point mass at the origin and let Ita,b) in) denote the in- 
dicator function of the real interval {a,b). Then, as 
t — > 00, the normalized position distribution ft{y) asso- 
ciated with ^Xt converges, in the sense of a weak limit, 
to the density function 



f{y) = coo6oiy) 



I, 



(-l/\/2,l/%/2) 



(y) 



7r(l- 2/2)^1 _22;2^ 



Y,c,y'- (18) 



In the above formula, the coefficients cqo, cq, ci and C2 
are given by 



Coo 



V2 1 



V2)i\ 



a2 



m\ 



— Re[(2 — v 2)(a2Q!4 -l- a^a^ — aia2 — cticts) 



-(3v2 — 4)aia4 — \2a20f^ 



Co 

Cl 



- + Re{a2a3 - a\aA\ 

|aip-|a4p 

-f Re(aia2 + Oi\az + Q;2cE4 -t- 01304) 

C2 = ^(|«lP + |«4p-|02p-|a3p) 

-l-Re[3aicJ4 -I- 0102 + 0103 
—0203 — dial — 030!)] . 



(19) 
(20) 



The derivation of Eg. ((18)) is due to the method by 
Grimmett et al. [Il|. 

For instance, if the 2cQRW is launched from the BeU 

state |$+) = 75(100) + |11)), then we obtain 

/(.)HV2-Wy) + ^^^^^-^i^f^44- (21) 

It is no accident that the coefhcient cqo = v^ — 1 of 
5q (y) in Eq. (|2T|) coincides with the value of the sum in 
Eq. (|12p . In this particular instance, we have 



Coo 



= ^p(x) = V2-l. 



xez 



A detailed study of the limit distribution of jXt for 
2cQRW will be presented in a forthcoming publication. 
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r ■^ fjh 

{e— ^■^(lS-(fc),M'e^(fc))T/,(fc)}- 

i=2 
{V,{27T),^-i2n))V,{2n)]+0{x-'). (A3) 

Therefore, 

i=2 
{V,{2tt), ^^{2TT))Vj{2TT)]r + o(x-^). (A4) 

At the other extreme, setting a; = 0, we obtain: 



p(0) 



2-n 3 



{Y,{V,{k),^l%k))V,{k)} 



J=2 



dk 

2^' 



(A5) 



APPENDIX B: PROOF OF THEOREM 2 



Proof. Let 



APPENDIX A: PROOF OF THEOREM 1 

Proof. Let ?ti > denote the highest order of the 
stationary points of <p{k). By the method of stationary 
phase [3|, Il8j applied to the first term and the fourth term 
in Eq. ([8]), we have 



271- ^ JU 

{e-"'=e'*^('=\yi(A:),*g"(fc))Vl(A:)} 



0{r 



27r 

1 , 

+1) as t ^ cxd; 



(Al) 



2^ ___ ji, 

{e~"'=e**('^-^('=»(V4(fc),*8'(fc))"^'4(fc)}7^ 



0{t ^5^+T) ast- 



2tt 

00. (A2) 



Since the eigenvalues corresponding to I2, h in Eq. 
([5]) are independent of fc, the values of I2 and J3 do not 
vanish as i — > c». Hence, I2 and I3 provide the only non- 
negligible contributions to the limiting amplitude. Thus 



p{x) 



r{e-"'=^(F,(fc),$p(fc))F,(fc)}^ 



By Lemma 1, as x — > 00, we obtain 



and let 



so that 



w{k) ^Y.^v,{ki^^mvAk). 

J=2 



c.= / e-"^W(fc)^, 
/o ^'^ 



VF(fc)=^c,e"'=. 



Since 



(W^(fc), W{k)) = Y. 11^-11' + E c.5^e^^^"'^'> 



x& 



and, by Theorem 1, 

Ep(^) = Eiic^ 



X\\ I 



we conclude that 

r2-ix 



f \w{k),w{k))^ = Y^wc^r 
= E^(^)- 



x& 



8 



On the other hand, 



(i¥(fc),M/(fc)) = ^i(y,(fc),$a^(fc))i^ 

J=2 



Therefore 



xGZ 



2t 3 



Ep(^)= / Ek^^W'*^wi^ 



i=2 



, dk 

2^' 



APPENDIX C: PROOF OF THEOREM 3 

Proof, (i) Assuming that 1 is the highest order of 
any stationary point of (p{k), we infer, by the method of 
stationary phase [3, [lg| , that the first term Ji and the 
fourth term I4 in Eq.® decay as follows: 



27r 



dk 



Ztt 
^ 0{t-i) ast->oo; (CI) 



2v 



{e-'''''e'*^'^^-'^^''^\V4ik),%%k))V4ik)} 



dk 



o{t- 



27r 
as t -> 00. (C2) 



Hence 



pM 



n2-n 3 



r Y^^v^(k),^i^[k))v,{k)^^f +o{t-'^) 



J=2 

p(0) + O(t-5 



(C3) 



(ii) If a; e \t^p'{kQ) — S,tip'{ko) +S], say, x = tip'{ko) + 
c, where \c\ < S, then the amplitudes of Eq.® can be 
expressed as follows: 



|git[(p(fc)-(p'(fco)fc]g-icfe/y^ ^^ec\Y^ 


dk 



_|_gj:t[26»-v(fc)-ip'(fco)fc]g-icfe/^ ^^\p^| 



27r' 



(C4) 



Since fco is a zero of f"{k) with multiplicity 1, fcg is 
a stationary point of </?(A:) — ip'{ko)k, 9 — kip'{ko) and 
26 — f{k) — ip'{ko)k of order 2, 0, and 0, respectively. 

By the method of stationary phase, as in [3|, [l8|, we 
infer that, as i — > 00, the first term /i in Eq. (|C4|) 
is the order of 0{t^^), while the other three terms in 



Eq. (|C4p are of order 0{t~^). Therefore, as t — > 00, 
tpt{x,j) = 0{t~i) where j = 1,2,3,4. 

Similarly, when x e [—tLp'{kf)) — 5, —tLp'{kf))+S\^ we get 
ipt{x,j) ^ 0(i^3)where j — 1,2,3,4. Hence, in either 
case, pt{x) = 0{t~3). 

(iii) For convenience, we may assume ip'{ko) > 0, then 
M — ip'iko). It suffices to treat the case when i > a; > 
((y5'(fco) + e)t, as other case can be treated similarly. The 
amplitudes of in Eq. ([5]) can be expressed as follows: 



{Mx,i),Mx,2),Mx,s),MxA))^ 



{e^*['^W-7'=l(-l/i,^)-l/i 



+e**^e-"'=(F2,*80"^2 



+e'''e~ 



^'(^^3,*80^3 

dk 



27r 



(C5) 



Note that ji:{(p{k) - f fc) < ^'{k) - ^'(fco) - e < -e, 
and ^(261 - ip{k) - fk) < -^'{k) - Lp' {k^) - e < -e. 
Thus, neither ip{k) — jk nor 26 — ip{k) — jk in. the ffist 
and fourth terms of Eq. (jCSp possess stationary point (s) 
of non-zero order. By the method of stationary phase 
in 0, [3, Ij ~ 0{t-^) as i -> 00, where j = 1,4. For 
the second and third terms in Eq. (jCSp . by appealing to 
Lemma 1, we have iptix, j) ^ 0{x^^) as t ^ cxd (therefore 
as a: — *■ 00), so Ij ~ 0{t^^) where j = 2,3. Hence 
p^{x)=0{t-^). 

(iv) It suffices to treat the case when \/t < x < 
{ip'{ko) — e)t, as the other case can be treated similarly. 
Note that -^ifik) - f A:) = tp'ik) - f . since -M + e < 
J < M — e, we see that ip{k) — |fc and 26 — ip{k) — jk 
both possess a stationary point of order 1. Thus, by the 
method of stationary phase in [1^, Ij ~ 0(i~2)^ where 



j = 1,4. Meanwhile, the other two terms in Eq. (|C5[) 
satisfy Ij ~ 0{x~^) where j = 2,3. Note that \/t < x < 
((/3'(fco) — e)^, so either Ij ^^ 0(i~2) or Ij ^ o{t^2) for 
j = 2, 3. Hence pt{x) = 0{t-^). 

(v) As in the proof of (iv), we have Ij ^ 0{t^^) for 
j = 1,4. Note that Ij ~ 0(1) as \x\ ~ for j = 2,3. 
Therefore 



Pt{x) 



(vi) As in the proof of (iv), it can be shown that Ij ^ 
0{t~i) for j = 1,4. Note that /j - Oi\x\-^) - 0(i"5) 
as |a;| approaches t^ for j — 2, 3, by Lemma 1. Therefore, 
pt{x)^0{t-^). 



f e-^^^)_^{V,{k),^l^{k))V,{k)^ 


2 

+ 0(t-5) 


Pix) + 0{t~i). 


(C6) 
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